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Electron transport through a quantum dot or single molecule coupled to a quantum oscillator is studied by 
the Keldysh nonequilibrium Green's function formalism to obtain insight into the quantum dynamics of the 
electronic and oscillator degrees of freedom. We tune the electronic level of the quantum dot by a gate 
voltage, where the leads are kept at zero temperature. Due to the nonequilibrium distribution of the 
electrons in the quantum dot, the spectral function becomes a function of the gate voltage. Novel spectral 
features are identified for the ground and excited states of nanomechanical oscillators that can be used to 
enhance the measurement sensitivity. 



I n the last decade, deeper understanding of a large class of systems important for electronic transport, such as 
I molecular systems, has been accomplished by the single resonant level model. Electron transport through 
I molecular single electron transistors has been a particularly hot topic in modern mesoscopic physics' 4 . The 
coupling of such a system to mechanical vibrations leads to important new physical phenomena, which have been 
demonstrated experimentally in many studies for molecular quantum dots, see Refs. 5-7, for example. In these 
experiments, the electron-phonon interaction is very strong and thus cannot be covered by perturbation theory. 
However, the quantum signatures of molecular vibrations have attracted much interest due to fundamental 
scientific challenges 8-10 and potential applications in nanoelectromechanical systems (NEMS) 11,12 . 

NEMS have become a promising prototype for measuring the quantum behavior of nanoscale objects, the- 
oretically as well as experimentally 12,13 . The main motivation to study NEMS is the large number of potential 
applications in nanoelectronic devices, such as the measurement of very small masses, weak forces, and ultra small 
displacements of atomic positions 111417 . NEMS are also important from a fundamental point of view for under- 
standing the transition between the classical and quantum regimes. This includes the observation, measurement, 
and control of quantum states of mesoscopic mechanical continuous variable systems, such as the quantum 
harmonic oscillator. Metallic leads coupled to the quantum dot have been successfully used to characterize the 
state of the mechanical resonator classically 18 and quantum mechanically 19 . When the device is scaled down to the 
nanometer regime, the resonance frequencies increase and at the same time the effective masses decrease. Related 
applications are discussed in a very recent review in Ref. 20. Various electromechanical properties of the devices, 
such as current, conductance, and shot noise, are proportional to the spectral function. Hence, the spectral 
features 2125 of the molecular systems are of great interest. In the latter studies, the spectral features have been 
discussed using the equation of motion and perturbative approximations for weak electron-phonon coupling. On 
the other hand, the spectral features of NEMS have not been discussed so far. 

In view of the above, we study in this work the spectral function of a single quantum dot connected to two 
metallic leads and a quantum harmonic oscillator. An arbitrary voltage is applied to tune the single electronic level 
of the dot. The quantum oscillator is in the ground state and the leads are assumed to be in equilibrium at zero 
temperature. We address very small to very large chemical potential differences between the Fermi levels of the 
left and right leads, where the response to both Fermi levels is included in the calculations. The position of the 
resonant level with respect to the chemical potential in the leads is thus affected by the displacement of 
the nanomechanical oscillator, which in turn affects the transport properties of the tunnel junction. Electronic 
transport properties of similar systems have been discussed for a single electron transistor 18 and a superconduct- 
ing single electron transistor 26 , both coupled to an oscillator. Both these studies have employed the quantum 
master equation approach within the weak coupling limit 20 . It has been observed that, even though the leads are at 
zero temperature, the oscillator behaves as if it would be coupled to a bath/reservoir of finite temperature related 
to the chemical potential difference. 
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Here we employ the nonperturbative Green's function approa- 
ch 1924 , as the recent experiments go beyond perturbation theory 
and thus a nonperturbative approach to the quantum master equa- 
tion (equation of motion) is required. The Green's function approach 
allows us to consider strong electron-phonon coupling and has been 
applied successfully for very small up to very large applied bias. In 
contrast, the equation of motion approach 18,25,26 treats both couplings 
as perturbations. In addition, this approach is valid only for high 
applied bias. We assume that quantum transport through the tunnel 
junction can be described by sequential tunneling. We note that the 
electron-oscillator interactions have been studied by scanning tun- 
neling microscopy 5,6,27 29 and spectroscopy 7,29,30 . These methods turn 
out to be more valuable tools than the inelastic tunneling spectro- 
scopy 31 33 . In Refs. 34-39 it has been shown that the position of the 
single quantum dot energy level and the tunneling width, T, strongly 
affect the Kondo spectrum, where a third terminal has been used as a 
tip detector to measure the Kondo peaks. Spectral features of the 
Kondo effect have also been observed in the presence of strong elec- 
tron-phonon interaction 40 42 . In the present work, therefore, we will 
consider whether the intriguing features of NEMS can be explored by 
the spectral function. 

Model formulation and the spectral function 

We consider a single quantum dot connected to two identical metal- 
lic leads and apply a gate voltage. The single mode of the quantum 
oscillator is strongly coupled to the electrons on the dot. We neglect 
the spin degree of freedom, the electron-electron interaction, the 
effects of finite electron temperature of the leads and the damping 
of the oscillator. The Hamiltonian H = H. 



dot-ph 



Hi 



leads < ^t]=L,R 
19.21-24 



H^-dot to be studied in the following is thus given by 

H dot . ph = [E 0 + ).{b f + b)]clc 0 + ntfb, (1) 

where E 0 is the energy of the single electron level on the dot and Cq , 
c 0 are the corresponding creation and annihilation operators. 
Moreover, X is the coupling strength between the electrons on the 
dot and an oscillator of frequency Q (setting H = 1 for simplicity) and 
y, b are the raising and lowering operators of the phonons. The 
remaining elements of the Hamiltonian are 



r]=L,Rj = l 



1 n t \ 
H ''~ dot ~ Vn +-( V '' v ' h ' c ° + c ° c ' hi ) ' 



(2) 



(3) 



where N is the total number of states in the lead and j counts the 
channels. We use the usual harmonic oscillator eigenfunctions 
¥„(*/€) when the electron is in the lead and *?„((* — x 0 )/€) when 
it is on the dot, where € is the zero point amplitude of the oscillator 
and x 0 = M/Q is the displacement of the oscillator equilibrium 
position due to the coupling to the electrons. The Hamiltonian thus 
takes the form 



H'- 



L 



^dot-ph 

vt 



(4) 



where H', 



do,- P h-L£o + (»+l/2)fi-A]<5„ y , H' ir [E, hJ +(n+l/2)n] 
P/2fl is an energy shift. 
In order to transform between the representations for the occu- 
pied and unoccupied dot we require the matrix elements 



difd n> rt, and A 



M. 



■x 0 )/£)V n (x/l) dx 



which may be simplified 43 as 



2 " 



-exp 



\n — m\ 



L . , I -x' 

min|m,n| I 2 



(6) 



where x = x 0 /£ = A/O. In addition, LJJ, (x) are the associated Legendre 
polynomials. 

The self-energy represents the contribution to the dot energy due 
to interaction between the dot and the leads it is coupled to. In the 
current model, with a single electronic state on the dot, there is no 
self-energy due to interaction. We use the wide-band approximation, 
where the self-energy contribution of each lead becomes energy 
independent and can be calculated as 



where 



Np^dE,, 



(7) 



(8) 



p n is the density of states in lead j; (L or _R), and g r nn (a>) is the 
uncoupled retarded Green's function in the leads (recall that k = 
1). The retarded self-energy can hence be written as 



(9) 



Here T, ; = 4rc | V n \ 2 p n and Y n is the damping factor {Y L = Y R = 
T). The retarded self-energy is now independent of the oscillator 
index. Similarly 



(10) 



We solve Dyson's equation using JJ^-dot as a perturbation. In the 
presence of the oscillator, the retarded and advanced Green's func- 
tions on the dot, with the phonon states in the representation of the 
unoccupied dot, may be written as 

(11) 



where g^ a (co) is the retarded/advanced Green's function on the 
occupied dot obtained as 



m) = \ co—E Q - 



m+i)Q + A + zT 



(12) 



The lesser Green's function in the presence of the nanomechanical 
oscillator including the dot and the leads is written as 

(13) 



Here £"*(«) is the lesser self-energy, given by 

S < (co) = S<(c 0 ) + S<(cq), (14) 

where the off-diagonal elements of the matrix (a>) are zero and 
the diagonal elements may be written as 



£ « < :«,,,( w ) = ! '') JdE n f„(E ri )B n 8(w-E ri -(n+^Q 



(15) 



Here e = (0— (n + ^JQ and/,,^) denotes the Fermi distribution 

functions of the left and right leads, which have different chemical 
potentials under the bias, and B n is the Boltzmann factor of the 
oscillator state. The index n determines the statistical occupation 
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probability of the phonon state | n) at finite temperature. Therefore, 
the accessibility of particular conduction channels is determined by a 
weight factor of the Boltzmann distribution function. For the present 
case of zero temperature the lesser self-energy may be recast in terms 
of the Heaviside step function 0(x) as 



■v +- Q- 



(16) 



where Ep is the Fermi energy on lead t] and the Kronecker delta, 8 n >0 , 
signifies that the nanomechanical oscillator is initially in its ground 
state, n = 0. All physical properties may be expressed in terms of the 
spectral function A(co), which is at the heart of this paper. The spec- 
tral function is related to the imaginary part of the retarded or 
advanced Green's function 2 1,22,25 



A(co)= -Im[G f ((»)] = 



-ImG fl ( 



which we obtain as 



2TM n>m M* Hm 



co-E 0 - 



(m+i)n+A]'+r 



■2' 



(17) 



(18) 



where n = n' = 0 for the ground state and 1 for the first excited state 
of the nanomechanical oscillator. The spectral function is then used 
to calculate the differential conductance via the relation 



1 = 



eT 
4tz 



dcoA(co)\f L (co) —fn{co) 



(19) 



where we have fixed the dot energy to E 0 = 0.5 and set the Fermi 
energy of the right lead to zero. The spectral function is proportional 
to the density of electron plus phonon states by the relation n(co) = 
A(to)/n, which gives information about the excitations (electron or 
hole) of the system. The spectral function obtained in Eqs. (17) to 
(19) is the main result of the present work and the starting point for 
the subsequent discussion. The form of the spectral function 
obtained here makes the electronic and oscillator properties apparent 
as it is the essential ingredient for theoretical considerations on pro- 
blems such as the high frequency behavior of nanomechanical oscil- 
lators, the current-voltage characteristics, quantum conductance, 
and shot noise phenomena 21-25,34-42 . 

Discussion 

Analytical examination of Eq. (17) demonstrates a resonant peak for 
the purely electronic degree of freedom plus any number, m, of 
phonons. This discrete structure confirms the quantum mechanical 




Figure 2 | Three-dimensional plot of the spectral function, as a function 
of oj and A, for the ground state (n = 0) of the nanomechanical oscillator. 

The same parameters as in Fig. 1 are used. 

aspect of the calculation 18,26 . We expect the side peaks to be sup- 
pressed in the absence of electron-phonon coupling through M m „ 
= 0 for m n. In this regime, the resonant peak is shifted in energy 
by a constant shift, A, as shown in Eq. (17). In all the subsequent 
figures the following parameters are used with the dimension of 
energy: Q, = 0.3, E 0 = 0.5, A = l.lfi, and T = 0.2fi. 

We discuss the spectral function in Fig. 1 for the ground state (n = 
0, dotted line) and the first excited state (n = 1, solid line) of the 
nanomechanical oscillator. We find n side peaks and a constant shift 
A consistent with Eqs. (4) and (17). This is the regime of strong 
coupling, A = 1.1Q, whereas in the limit of weak coupling, A = 
0.5fi, we find a single side peak for the excited state. The amplitude 
of the side peaks becomes larger and the peaks sharper when the 
oscillator starts in the excited state instead of the ground state. We 
find a nontrivial feature for the first excited state of the nanomecha- 
nical oscillator, see the solid lines in Fig. 1: The first side peak is 
suppressed with a very small amplitude as compared to all other side 
peaks for strong electron-phonon coupling, which is not the case 
when the nanomechanical oscillator starts in its ground state (dotted 
lines in Fig. 1). 

In order to explore this nontrivial behavior, we show in Figs. 2 and 
3 three-dimensional plots of the spectral function, as a function of (U 
and the coupling strength A, for the ground (n = 0) and excited (n = 
1) states, respectively. We observe that valleys appear when the values 
of the excited state and integer values coming from the probability 
density function in Eq. (17) are the same. Closer examination of this 
equation shows that for m = 1 the side peaks are suppressed for 
different coupling strengths. Such a behavior is not possible in clas- 
sical dynamics 18,26 or any other perturbative approach used for 
molecular systems 21-25 . Finally, we show the probability density 
P(x) =M m ^„M* calculated by Eq. (6) in Fig. 4 as a function of the 




Figure 1 | Spectral function of the ground (n = 0, dotted line) and first 
excited (n = 1, solid line) states of the nanomechanical oscillator. The 

employed parameters (dimension of energy) are: Q = 0.3, E 0 = 0.5, A = 
l.lfi and T = 0.2CX 



Figure 3 | Three-dimensional plot of the spectral function, as a function 
of co and A, for the first excited state (n = 1) of the nanomechanical 
oscillator. The same parameters as in Fig. 1 are used. 
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2 3 
x (Dimensionless coupling) 

Figure 4 | Probability density as a function of the dimensionless coupling 
strength x for the ground (n = 0, dotted line) and first excited (n = 1, 
solid line) states of the nanomechanical oscillator. The same parameters as 
in Fig. 1 are used. 

dimensionless coupling strength x to confirm the suppression of the 
side peaks. Results for the ground state and the first excited state of 
the nanomechanical oscillator are given by dotted and solid lines, 
respectively. From the solid line it is clear that when the oscillator 
starts in an excited state the curve approaches zero for n values of the 
coupling. Using Eq. (6), we obtain for the probability 



(20) 



P(*)x[L?(p 
which can be simplified using for x > 0 the identity 4 ' 

!"(*) = -L e x/2 x -(2« + l)/4„(2»-l)/4 C0 J 2 y^_ 

y/ll V 

+o(«' 2m - 3 »/ 4 ) 



(2m+ 1)71 



(21) 



to yield 



P(x)oz 



e x ^ 4 {x 2 /2)- 1/4 cos 



y/7x-l 



(22) 



No such behavior has been reported for classical and perturbation 
treatments. The reason is that in the classical mechanics and per- 
turbation approaches the oscillator does not couple strongly enough 
to the electrons on the dot. The nontrivial behavior only appears 
because a nonperturbative approach has been used in the present 
work. As compared to the existing literature, our nonperturbative 
approach is valid for small to large applied bias and can be employed 
for studying various problems in molecular systems 1 - 4,21-25 , 
NEMS 8 - 14 , and shuttling transport 45 - 52 . 

The model considered in this work represents a nanomechanical 
oscillator coupled to a resonant tunnel junction, which can be 
directly used as small mass detector, sensor of ultra weak forces, 
and particularly as displacement detector 8 " 12,19,20 . We propose to 
move from the quantum ground state of the nanomechanical oscil- 
lator to an excited state, because the resonant electronic peak will be 
suppressed and the mechanical side peaks will become sharper with 
larger amplitude. Due to the enhanced amplitude for strong electron- 
oscillator coupling, the quality factor is enhanced, which increases 
the measurement sensitivity in the quantum regime 53-58 . The 
required excited state of the nanomechanical oscillator can be 
achieved by controlling the temperature as well as by applying an 
external gate. The application of NEMS often suffers from difficulties 
to control the damping in the systems. Various ways to overcome 
damping losses and, therefore, to improve the quality factor have 
been put forward both theoretically and experimentally, such as 



the parametric resonance 53,54 , mechanical signal amplifying 55 , 
sideband cooling 56,57 , and the quantum behavior of macroscopic 
systems 58 . The enhanced amplitude of the sidebands in the presence 
of strong electron-oscillator coupling for the excited state, as demon- 
strated in the present work, gives rise to an alternative and very 
efficient approach for increasing the measurement sensitivity of 
NEMS. 

We have analyzed the spectral features of a quantum dot coupled 
to a nanomechanical oscillator by using the nonequilibrium Green's 
function approach, i.e., without treating the electron-phonon coup- 
ling as a perturbation. We have derived analytical expressions for the 
self-energies of the dot and the leads as well as the full Green's 
function. This enables us to investigate the effects of the coupling 
of the electrons to the oscillator on the dot and the tunneling rate of 
the electrons. Our numerical results show additional peaks in the 
spectrum for strong coupling to the single phonon mode. We find 
nontrivial dynamics of the nanomechanical oscillator, as reflected by 
a suppression of the first side peak in the spectrum when the oscil- 
lator starts in its first excited state. This effect could not be revealed by 
a classical or perturbation treatment of the oscillator. While both for 
the ground and excited states the side peaks are well resolved, they are 
enhanced in amplitude and sharper for the first excited state with 
respect to the ground state. Our method goes beyond perturbation 
theory and thus can be applied from small to large applied bias. It 
should turn out to be fruitful for various problems in molecular 
systems, NEMS, and shuttling transport phenomena. The theoretical 
framework developed in the present work for describing the ground 
and excited states of NEMS provides both reliable predictions of 
experiments and improved understanding of the circumstances 
under which nanomechanical systems are governed by quantum 
behavior. 
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